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Abstract. Wc study the L- functions associated to Siegel modular forms (equivalently, 
automorphic representations of GSp(4,Aq)) both theoretically and numerically. For the 
L-functions of degrees 10, 14, and 16 we perform representation theoretic calculations to 
cast the Langlands L-function in classical terms. We develop a precise notion of what it 
means to test a conjectured functional equation for an L-function, and we apply this to the 
degree 10 adjoint L-function associated to a Siegel modular form. 



1. Introduction 

L-functions are special functions that arise in representation theory and in several areas 
of number theory. From the viewpoint of analytic number theory, L-functions are Dirichlet 
series with a functional equation and an Euler product. From the point of view of represen- 
tation theory, L-functions arise from automorphic representations of a reductive group over 
the adeles of a number field. 

The two points of view offer distinct benefits. Representation theory, via the Langlands 
program [TU], provides a framework for understanding how L-functions arise, as well as the 
connections between various mathematical objects. L-functions considered as objects of an- 
alytic number theory are suitable for concrete exploration and testing of conjectures, for 
example they can be evaluated on a computer to check the Riemann hypothesis. Unfortu- 
nately, it can be quite difficult to translate Langlands L-functions into this setting, which 
limits the ability to do explicit calculations and test conjectures. In this paper we make such 
a translation and perform computer calculations with the results: for a particular Siegel 
modular form F, we calculate factors Lp{s,F,p) and ep{s,F,p) {p < oo) for six choices of 
p (dimensions 4, 5, 10,14, 16) using the Langlands parameterization of the discrete series. 
Using these calculations we provide numerical evidence that the L-function of degree 10 
satisfies a functional equation. 

The L-functions we consider here arise (in the classical setting) from holomorphic Siegel 
modular forms on Sp(4, Z), see Section [2j The same L-functions arise from automorphic 
representations of PGSp(4, A), see Section [3] (as in |2], at the archimedean place such an 
automorphic representation is a holomorphic discrete series representation and at the nonar- 
chimedean places it is a spherical principal series representation). The Langlands program 
predicts the existence of an infinite list of L-functions associated to a Siegel modular form. 
In our particular case the first two L-functions are known as the spinor and the standard 
L-function, and have degree 4 and 5, respectively. Due to Andrianov pQ, Shimura [2T], 
Bocherer [3] and others, these L-functions are fairly well-understood; for Siegel modular 
forms on Sp(4, Z), they are known to be entire functions that satisfy a functional equation. 
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The next case is the adjoint L-function, which has degree 10. It has not been shown that 
this L-function is entire and satisfies a functional equation. (The automorphic represen- 
tations associated to holomorphic Siegel modular forms are not generic, so a technique as 
in [11] is not applicable.) Providing evidence for the conjectured functional equation, via a 
computer calculation, is one of the goals of this paper. This is made precise in Theorem 14. 2[ 
which gives a test for the functional equation and quantifies the probability that the test 
will yield a false positive. 

A substantial part of this paper is a translation from the perspective of representation 
theory to the viewpoint of analytic number theory. Selberg [T^] gave a set of axioms for 
what is now called the "Selberg class" of L-functions. We will call L-functions in this class 
"Selberg L-functions," and are to be compared with Langlands L-functions - those that arise 
from automorphic representations. It is a standard conjecture that all Selberg L-functions are 
Langlands L-functions and that all primitive Langlands L-functions are Selberg L-functions. 
In this paper we translate Langlands L-functions into Selberg L-functions. We now describe 
Selberg L-functions in more detail. 

A Selberg L-function L(s) is given by a Dirichlet series L{s) = ^ where oi = 1 and 

the series converges in some half-plane. We assume a Ramanujan hound on the coefficients: 
On = 0{7f) for any e > 0. Moreover, it has a meromorphic continuation to the whole 
complex plane with at most finitely many poles, all of which are on the line 3f?(s) = 1. L{s) 
can be written as an Euler product L{s) = Ylp-^piP'^)'^ where the product is over the 
primes, and Lp is a polynomial with Lp{0) = 1. Additionally, there exist Q > 0, positive real 
numbers Ki,. . . and complex numbers with non-negative real part fii,. . . ,fin, such that 
A(s) := Yl r('^j'5 + /ij) ■ L{s) is meromorphic with poles only arising from the poles of L{s) 
and satisfies the functional equation A(s) = £:A(1 — s) where = 1. The number d = 2^ kj 
is the degree of the L-function. An alternate way of thinking about the Ramanujan bound 
is that I Op I < d. 

Indeed, for the L-functions considered here there do not exist results in the literature which 
would allow a non-expert to translate the L-function data from representation theoretic 
language into a form involving a Dirichlet series. Thus, we give a brief introduction to the 
aspects of the Langlands program that are relevant to our calculations, and describe how to 
translate the Langlands L-functions we consider here into Selberg L-functions. The results 
of those calculations are summarized in Proposition 12.11 

In Section [2] we describe the Siegel modular forms which give rise to the L-functions 
considered here, and we describe these L-functions in the classical language. In Section [3] 
we describe Langlands L-functions and how to translate the degree 10 L-function considered 
here into a form which can be evaluated on a computer. In Section H] we provide evidence 
for the conjectured functional equations and also briefiy address the problem of accurately 
evaluating L-functions for which only a few of the local factors in the Euler product are 
known. We also provide criteria to measure the strength of the evidence. 

We would like to thank Sharon Garthwaite for her helpful comments on the paper and the 
referee for helping improve the paper. 

2. Siegel modular forms and their L-functions 

We recall the definition and main properties of Siegel modular forms on Sp(2n, Z) and we 
describe the two simplest L-functions associated to them. 
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2.1. Siegel modular forms. Let 0^ be the zero matrix, En '■= .-' j and J„ : = 

( \ ) . Denote the group of symplectic simihtudes by GSp'''(2n, M) := {a G GL(2n,M) : 
*q;J„q; = r{a)Jn,r{a) G M, r(a) > 0)}, where r{a) is called the similitude of a. Define the 
Siegel modular group of genus n by Sp(2?7,, Z) := {7 G GSp^(2r?,, M) fl M(2n, Z) : r(7) = 1}. 
Let W := {Z = X + lY : X,Y e M(n,R),*Z = Z,Y > 0} denote the Siegel upper half 
space, that is, symmetric matrices in M{n, C) with positive definite imaginary part. 

Recall, a holomorphic function F : "H" — )■ C is a Siegel modular form of genus n and 
weight k if for all a = ( ^ ^ ) G Sp(2n, Z) it satisfies the transformation property 

F{Z) = {F\ka){Z) := r(a)"'=-"("+i)/2 det(CZ + D)-''F{{AZ + B){CZ + D)-^). 

If 77- = 1 then F must satisfy an additional growth condition. 

We shall denote the space of weight k genus 2 Siegel modular forms by Mfc(Sp(4, Z)). 
In genus 2, we can express the expansion of a Siegel cusp form as 

F{Z)= Yl aF{n,r,m)q\^r 

r,n,m£Z 
r^— 4mn.<0 
n,m>0 

where [n, r, m] is the positive definite binary quadratic form nX^ +rXY +mY'^ of discriminant 
- 4mn and q = e^'^" {z G H^), q' = e^'^''^ (w G H^), and C = e^''*^ (r G C). In particular, 
we are examining L-functions associated to modular forms not in the Maass space, i.e., 
not in the image of the Saito Kurokawa lift, i.e., whose L-functions are primitive (at least 
conjecturally). The first such form occurs in weight = 20 and is computed in [22] . 

There is a theory of Hecke operators acting on the space of Siegel modular forms; we 
denote the nth Hecke operator by T{n). The Hecke eigenvalues for a Siegel modular form can 
be computed explicitly from its Fourier coefficients, but this is computationally expensive. 
Let F be a Hecke eigenform; i.e., suppose that for each n there exists a Xpin) so that 
F\T{n) = Air(n)F. For the weight 20 Siegel cusp form F that is not a Saito-Kurokawa lift, 
we will use Xf{p) and Xpip^) for p < 79 in our experiments. These data are computed in [15]. 

2.2. L-functions associated to Siegel modular forms. There are two well-known L- 
functions attached to Siegel modular forms on Sp(4, Z), called the spinor L-function and the 
standard L-function. These have been studied by Andrianov [T], Shimura [21], Bocherer [3] 
and others. Formulas for those L-functions in the genus 2 case are given in Proposition 12.11 

To each genus 2 eigenform, F, one can associate, for each prime p, a triple (ao,p, ai,p, a2,p) 
of nonzero complex numbers - it is in these terms that our L-functions are expressed. The 
entries of the triple are called the Satake parameters of F. 

In genus 2 it is rather straightforward to compute the Satake parameters of a form, given 
the Hecke eigenvalues of the form. By using an explicit description of Satake isomorphism as 
found, for example, in [18], write the Euler factor of the spinor L-function as a polynomial 
whose coefficients are in terms of Xf{p) and Xp^p"^). To compute the Satake parameters, one 
finds the roots of this polynomial. 

We rescale the Satake parameters to have the normalization \aj\ = 1, aQaia2 = 1, which 
is possible since the Ramanujan bound is a theorem for Siegel modular forms [24]. This 
corresponds to a simple change of variables in the L-functions, so that all our L-functions 
satisfy a functional equation in the standard form s o 1 — s. 
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In Section [3] we describe the procedure for determining the L-functions associated to an 
automorphic representation, and give reasonably complete details for the spinor, standard, 
and adjoint L-functions of genus 2 Siegel modular forms. The results of those calculations 
are summarized in the following proposition. For completeness we also include the results 
of similar computations carried out for two more L-functions: ones associated to a specific 
14-dimensional representation and 16- dimensional representation. 

Proposition 2.1. Suppose F G Mfc(Sp(4, Z)) be a Hecke eigenform. Let ao,p, cxi^p, a2,p 
he the Satake parameters of F for the prime p, where we suppress the dependence on p in 
the formulas below. For p G {spin, stan, adj, pi4, pig} we have the L-functions L{s,F,p) := 
rip p„me QpiP"^ where 



Qp{X, F, spin) 
Qp{X, F, Stan) 
Qp(X,F,adj) 



Qp{X,F,p 
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1^1 — aoX){l — aoaiX){l — aoa2X){l — aoaia2X), 
[1 - X)(l - aiX){l - a^^X){l - a2X){l - a^^X), 
[1 - X)2(l - aiX){l - ayX){l - a2X){l - a^^X) 

(1 — aia2X){l — ai'^a2X){l — ai«2^^X)(l — a^^a2^X), 
[1 - X)\l - aiX){l - a2X){l - ayX){l - a^^X) 



(2.1) 



(1 - alX){l - alX){l - a^^X){l - a^^X) 

(1 — aia2X)(l — aia2^X)(l — ai^a2X)(l — ^a^^'^X), and 

aoXY{l — aoaiX)^(l — aoa2^)^(l — «oai«2^)^ 



aoayX){l — aoa2^X){l — aoa^X)(l — aoa^X) 

aoaia2^X){l — aoa^^a2X){l — aoala2X){l — aoaialX). 

give the L-series of, respectively, the spinor, standard, adjoint, degree 14 and degree 16 L- 
functions. These L-functions satisfy the functional equations (the last three conjecturally 
satisfy the functional equations): 

A(s, F, spin) := (47r2)-T(s + |)r(s + k- |)L(s, F, spin) 

= (-l)'=A(l-s,F,spin), 

(4Tr^/Y'T{ls)r{s + k- 2)r(s + k- l)L(s, F, stan) 

A(l — s, F, stan). 



A(s, F, stan) 
A(s, F, adj) 



[167r^)-T(i(s + l))'r(s + l) 
X r(s + A; - 2)r(s + k- l)V{s + 2k - 3)L(s, F, adj) 
= A(l-s,F,adj), 
A(s, F, Pm) := (2V)"^r(s/2)2r(s + l)r(s + k- 2)T{s + k - 1) 

X r(s + 2k- 4)r(s + 2k~ 3)T{s + 2k - 2)L(s, F, p^) 
= A(l — s, F, P14), and 

A(s,F,pi6) := (2V)-^r(s + i)'r(s + k- |)r(s + k-ffT{s + k- i) 
X r(s + 2A;- |)r(s + 2/e- |)L(s,F,pi6) 
= -A(l-s,F,pi6). 



(2.2) 
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3. Langlands L-functions 

The Euler products in the previous section arise as Langlands L-functions attached to 
automorphic representations of GSp(4, A) generated by the Siegel modular form F, see 
[2]. In general, this procedure involves the local Langlands correspondence, which is now a 
theorem for GSp(4); see [9]. However, since we are only interested in full level Siegel modular 
forms, the mechanism simplifies considerably. We shall briefly describe how to obtain the 
local factors in the non-archimedean and the archimedean case. 

The non-archimedean factors. Let ao, O-i, 0L2 be the Satake parameters of F at a finite place 
p, normalized as above, so that Oi^a\a2 = 1. These determine a semisimple conjugacy class 
in the dual group Sp(4, C), represented by the diagonal matrix 

A^^ = Diag (ao, aotti, ao«2, «ottia2) • (3.1) 

(One has to carefully go through the definitions of the local Langlands correspondence to 
see this; see Sect. 2.3 and 2.4 of jUj.) Let p : Sp(4,C) — )■ GL(?t,,C) be a finite-dimensional 
representation of the dual group. The local L-factor attached to the data ao, ai, 012 and p is 
given by 

The three smallest non-trivial irreducible representations of Sp(4, C) are the four- dimensional 
"spin" representation (which is simply the inclusion of Sp(4,C) into GL(4,C)), the five- 
dimensional "standard" representation (described explicitly in appendix A. 7 of [16]), and 
the ten-dimensional adjoint representation adj on the Lie algebra sp(4,C). Calculations 
show that the resulting L-factors are given as follows, 

Lp(s, F, spin) = (5p(p"^ F, spin), 
Lp(s, F, stan) = Qp{j)~^ , F, stan), 
Lp(s, F, adj) = Qp{p~\ F, adj), 

with the factors Qp as in the previous section. There are also corresponding local e-factors 
ep{s,F,p), which for unramified representations are all constantly 1. 

The archimedean factors. The real Weil group IVr is given by Wm. = U jC^ with the 
rules = —1 and jcj~^ = c (see [23] (1.4.3)). The commutator subgroup is C , the 
set of complex numbers with absolute value 1. There is a reciprocity law isomorphism 

rM : M>< ^ Wl\ (3.3) 
-l^jS\ 

M>0 9 X I — > y^S\ 

Let I ■ I be the usual absolute value on M, and let || • || be the character of defined by the 
commutativity of the following diagram, 

(3.4) 
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(see [23] (1.4.5)). Hence, \\z\\ = \z\'^ for 2 G C^, where | ■ | denotes the usual absolute value on 
C. The character || ■ ||* is denoted by Us, for a complex number s (see [23] (2.2)). There are 
L- and e-factors attached to characters of (see [23] (3.1)). The correspondence between 
characters of and characters of M^, and the associated L- and e-factors, are given in the 
following table. 



char, of Wr 


char, of 


L-factor 


e-factor 




z\'^\ j ^ 1 




7r-(-+<)/2r(^) 


1 




J ^ -1 


X I—)- sgn(a;) X * 


^-(s+t+l)/2Y(^s+t+l-^ 


i 



(3.5) 



Besides one-dimensional representations, the only other irreducible representations of 
are two-dimensional and indexed by pairs {i, t), where £ is a positive integer and t G C. The 
representation attached to {i,t) is (pi^t, given by 



iff 

The associated L- and e-factors are 
Lis 



^2£ gi- 









1 



(3.6) 



(3.7) 



Now, to a Siegel modular form F of weight k there is attached the four-dimensional repre- 
sentation of Wm. given by 

9^(fc-i,fc-2) := Vi,o © <^2k-3,o (3.8) 

(this is really the parameter of a holomorphic discrete series representation with Harish- 
Chandra parameter [k — l, k—2), see [S]). The image of this parameter can be conjugated into 
the dual group Sp(4, C). Given a finite-dimensional representation p : Sp(4, C) — )■ GL(n, C), 
we compose p with the representation (13. 8p and obtain an n-dimensional representation of 
Wk. By [13], this representation can be decomposed into one- and two-dimensional irre- 
ducibles. The product of the L-factors (resp. e-factors) attached to these irreducibles is by 
definition Loo{s,F,p) (resp. eoo{s, F, p)). Calculations show that 

Loo(s,F,spin) = 4(27r)-2^-'=+ir(s 

Loo(s, F, Stan) = 2-2^-2'=+5 7^-^^/2-2^+3 r(is)r(s + k- l)r(s + k-2), 



l)T{s + k-l), 



Loo{s,F, adj) = 2 



-4s-3A;+9^-5s 



-3^-+4r(i(s + i))2r(s 



xT{s + k- 2)T{s + k- l)V{s + 2k- 3), 



and 



£:oo(s,-F, spin) = (-1) £00(5, stan) = 1, 5oo(s, adj) = 1. 

We see that, up to an irrelevant constant, the archimedean L-factors coincide with the F- 
factors in ([22]). 

The global L-function. Having defined all local factors, the global L-function attached to F 
and a finite-dimensional representation p : Sp(4, C) — GL(r;,, C) is given by 

A(s,F,p)= J]Lp(s,F,p). 

p<oo 
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Up to a constant, this definition coincides with the Euler products defined in (12.11) . By 
general conjectures, the global L-function, which is convergent in some right half-plane, 
should have meromorphic continuation to all of C and satisfy the functional equatioiu 

A(l-s,F,p)=£(s,F,p)A(s,F,p), where e{s, F, p) = H e^is, F, p). (3.9) 

p<oo 

Note that in our case e{s,F,p) = eoo{s, F, p). Hence, the functional equations (12. 2p are all 
special cases of the general conjectured functional equation (13. 9p . 

4. Checking the functional equation 

As mentioned in the introduction, the degree 10 adjoint L-function associated to a Siegel 
modular form has not been proven to satisfy a functional equation. We develop a method of 
checking a conjectured functional equation, and in Theorem 14.21 we provide a quantitative 
result that estimates the probability that this test could yield a false positive. 

The main idea behind our method of testing a functional equation is that an L-function 
can be evaluated, at a given point and to a particular accuracy, using finitely many of 
its Dirichlet series coefficients. That evaluation makes fundamental use of the functional 
equation. Furthermore, this can be done in more than one way. The consistency of those 
calculations provides a check on the functional equation. We quantify the "probability" that 
the calculations are accidentally consistent by viewing the coefficients of the L-function as a 
random variable. 

In the next section we describe the approximate functional equation and use it to evaluate 
an L-function. Then in Section Hl2l we elaborate on the ideas from [8] to develop our method 
to check the functional equation for the degree-10 Euler product associated to a Siegel 
modular form. 

4.1. Smoothed approximate functional equations. The material in this section is taken 
directly from Section 3.2 of [T7] . 
Let 



oo 



n=l 

be a Dirichlet series that converges absolutely in a half plane, 3?(s) > cti. 
Let 

A{s)=Q^ (^n^jS + X,)^ L{s), (4.2) 

with Q, Kj G M^, 3^(Aj) > 0, and assume that: 

(1) A(s) has a meromorphic continuation to all of C with simple poles at si, . . . , se and 
corresponding residues ri, . . . , r^. 

(2) A(s) = £:A(1 - s) for some e G C, \e\ = 1. 

(3) For any < o's, L{a + it) = O(expt^) for some A > 0, as \t\ — j- oo, < cr < as, 
with A and the constant in the 'Oh' notation depending on a2 and a^. 



^Tlie local e-factors also depend on the choice of a local additive character. We are assuming a standard 
choice and hence do not reflect it in the notation. The global e-factor is independent of the choice of global 
additive character. 
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To obtain a smoothed approximate functional equation with desirable properties, Rubin- 
stein [T7] introduces an auxiliary function. Let (7 : C — ?■ C be an entire function that, for 
fixed s, satisfies 

|A(^ + s)g{z + s)z~'^\ 

as \'^z\ — 7- 00, in vertical strips, —Xq < ^z < Xq. The smoothed approximate functional 
equation has the following form. 

Theorem 4.1. For s ^ {si, . . . , sg}, and L{s), g{s) as above, 

\^-^ s — si. ^-^ 

^k=l " n=l 

(4.3) 

n=l ^ 

1 rv+ioo °- 

f,{s,n):=— / l[T{K,{z + s) + Xj)z~'g{s + z){Q/nydz 

^'^'^ Ju-ioo 
1 pu+ioo ^ 

f,{l-s,n):=— l[T{K,iz+l-s)+T,)z-'g{s-z){Q/nydz (4.4) 



where 



with V > max {0, —^{Xi/ki + s), . . . , —^{Xa/Ka + s)}. 

We assume L{s) continues to an entire function, so the first sum in (14. 3 p does not appear. 
For fixed Q, k. A, e, and sequence bn, and g{s) as described in the next section, the right side 
of (14.31) can be evaluated to high precision. 

We illustrate the approximate functional equation with an example. Our examples use the 
genus 2 Siegel modular form which is the unique weight 20 eigenform, F, that is not a Saito- 
Kurokawa lift. We consider the degree 10 adjoint L-function associated to F, which we denote 
by Lp^iQ. Conjecturally it satisfies the functional equation given in (12. 2p in Proposition 12.11 

It is convenient to instead evaluate the "Hardy function" Zp^io associated to Lp^iQ. This 
function is defined by property that Z^io(| + it) is real if t is real, and |Z^io(| + ^^)| = 
|-^F,io(| + it)\- We use L and Z interchangeably in our discussions. 

If we let g{s) = 1 and s = 1/2 + i then (14.31) gives 

ZfMi + = 1-15426 + 0.778012 62 + 0.50246 63 + 0.33776 64 + 0.235813 65 

+ ■ ■ ■ + 0.0000142432 632 + 0.0000132692 fegs + ■ ■ • 

+ 2.8771 X 10"^ 6149 + 2.7402 x 10"^ 6150 + ■ ■ ■ • (4.5) 

Here and throughout this paper, decimal values are truncations of the true values. The 
numerical calculations were done in Mathematica 7.0.1.0 on a Dell Inspiron 9300 laptop 
running RedHat Linux. 

If instead we let g{s) = e~^**/^ and keep s = 1/2 + i then (14. 3 p gives 

ZfMI + ^) = 1-3044 + 0.678149 62 + 0.314111 63 + 0.12853 64 + 0.0341584 65 
+ ■ ■ ■ + 0.0000147237 632 + 0.0000123925 633 + ■ ■ • 

- 1.28515 X 10"^ 6149 - 1.22359 x 10"^ 6150 + • • • . (4.6) 
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1 

2 


{+1 











2.148764 
± 0.00016 


3.084662 
± 0.00046 


3.263120 
± 0.0044 


-0.403124 
± 0.071 


0.446949 
± 1.48 


i 

4 


2.148757 
± 0.000027 


3.084643 
± 0.00011 


3.262960 
± 0.0013 


-0.405569 
± 0.023 


0.396311 
± 0.50 


1 
2 


2.148743 
± 0.00021 


3.084617 
± 0.00034 


3.262768 
± 0.0019 


-0.407940 
± 0.018 


0.356202 
± 0.21 


3 
4 


2.148744 
± 0.00014 


3.084617 
± 0.00025 


3.262767 
± 0.014 


-0.407989 
± 0.0019 


0.355043 
± 0.16 


1 


2.148772 
± 0.00039 


3.0846503 
± 0.00037 


3.262906 
± 0.0011 


-0.406974 
± 0.0056 


0.365212 
± 0.039 


2 


2.148146 
± 0.0087 


3.084355 
± 0.0040 


3.262411 
± 0.0069 


-0.408305 
± 0.019 


0.361331 
± 0.071 


3 


2.296591 
± 2.55 


3.108788 
± 0.42 


3.277819 
± 0.25 


-0.391079 
± 0.26 


0.388505 
± 0.38 


Table 4.1. Values obtained for Zpiq{\ + iT, (3) using equation 


IB]) with test 



function g{s) = e , and using the known Satake parameters for F for p < 79. 



To obtain a numerical value for Zf,io{\ + i) we need to know the coefficients The 
calculations in [T3] provide the Satake parameters of F for the primes p < 79, so Proposi- 
tion [2]T] gives the local factors in the Euler product for p < 79. Expanding the product gives 
values for infinitely many 6„, including all n < 82, all composite 83 < n < 79^, etc. Using 
the known 6„ gives an approximation to Zi? io(| + and the numbers in (14. 5 p or f l4.6p make 
it seem plausible that the contribution of the unknown bn is small, but we wish to make this 
precise. This has two ingredients: a bound on bn and a bound on the terms which appear 
in the approximate functional equation. For the unknown 6„ we assume the Ramanujan 
bound, so for example \bp\ < 10 if p is prime. For the contribution from the terms in the 
approximate functional equation, we use the bound in Lemma 14.31 for large n, and directly 
calculate for smaller n. See Section 14.31 for more details. 

The results for (14. 5 p and (14. 6 p are, respectively 

ZfMI + = 3.084662 ± 0.00047 

ZpAoik + = 3.084649 ± 0.00056. (4.7) 

The values in (14.70 are consistent with each other. We view this as a confirmation of the 
conjectured functional equation for Lp^iQ. 

We summarize the results of similar calculations, for various s and functions g, in Table l4?n 
Each column of the table corresponds to a value for s, and each row corresponds to a 
function g{s) = e~*^* in Theorem 14. 1[ Scanning down each column shows that the values 
are consistent, which gives a check on the functional equation for Lf.io- We make this more 
precise in the next section. 

4.2. Numerically checking the functional equation. We wish to check that the adjoint 
L-function given by the Euler product (12.10 has an analytic continuation which satisfies the 
functional equation. This requires that we evaluate the function outside the region where the 
Euler product converges, and also check that these values are consistent with the functional 
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equation. In what follows we fix the genus 2 Siegel modular form to be the unique weight 
20 eigenform, F, that is not a Saito-Kurokawa lift. Recall that we can use the calculations 
in [15] to determine the Satake parameters of F for the primes p < 79. 

Let g{s) = g{s,(3) = e~^^^ in Theorem 14.11 This meets the conditions of the theorem 
if 1/3 1 < f Xl^j- We use equation f l4.3p to test the functional equation. This cannot be 
done in a naive way, because A(s), as given by the right side of f l4.3p . automatically satisfies 
A(s) = £:A(1 — s). Instead we exploit the fact that the right side (14. 3p . with our choice of (?, 
has /3 as a free parameter. 

Rewrite (14. 3p as A(s) = g{s)~^T{s, Q, k, A, e, {&„}, (3), and let 

L{s,l3) = Q-' (j[ri^,s + X,)^ gis,pr'Tis,Q,K,X,e,{br,},/3). (4.8) 

Our test for the functional equation of L{s) is that L{s,f3) is independent of /3. That is, we 
check the consistency equation 

Z{s,f3i)-Zis,f32) = 0. (4.9) 
For example, using (14.50 and (14. 6 p gives (14. 9 p in the form 

ZfMI + hO) - ZfMI + h^f^) 

= 0.150138 - 0.0998628 62 - 0.188349 63 - 0.20923 h - 0.201655 65 
+ ■ ■ ■ + 4.80503 X 10~^ ^82 - 8.76677 x 10^^ bss + ■ ■ ■ 
- 1.57286 X 10"^ 6149 - 1.49761 x 10"^ 6150 + ■ ■ • 

= 0. (4.10) 

As described immediately before equation (14.70 . we can estimate the contribution of the 
bn which are not known. The result is 

^f,io(i + ^,0)-Zp,io(| + ^,3/2) 

= 0.150138 - 0.0998628 62 - 0.188349 63 - 0.20923 64 - 0.201655 65 

+ ■ ■ ■ + 4.80503 X 10^^ 682 + 1.49761 x 10^^ 6150 + ■ ■ ■ 

= e X 0.00077, (4.11) 

where |6| < 1. 

Now we can explain our method of evaluating the strength of (14. lip as a test of the 
conjectured functional equation. We wish to quantify the intuitive notion that it is unlikely 
for (14. lip to be true just by chance, because the coefficients of the 6„ are large compared to 
the right side of the equation. We do this by considering the 6„ to be random variables, and 
furthermore we make some assumptions about their probability density functions. This, of 
course, requires some justification which we now provide. 

L-functions naturally fall into families fi2\ [7] , and the collection of L-functions in a family 
can be modeled statistically. For example, for the family of GL(2) L-functions, each co- 
efficient bj has a particular distribution. The distribution of bp, for p prime, tends to the 
Sato- Tate distribution as p — )■ 00, and bn, bm are uncorrelated if {n,m) = 1. See [201 E] for 
details. 

For other families, there are other distributions, see [13] for several examples. These dis- 
tributions are the distributions of traces of matrices in a compact group, weighted according 
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to Haar measure. For the Siegel modular forms we consider here, the Hecke eigenvalues 
are expected to be distributed according to an Sp{4:, Z) analogue of the GL{2) case. This 
leads to a conjecture for the distribution of the Dirichlet series coefficients of the degree-10 
L-function we are considering here. (See additional comments at the end of this section.) 

Thus, over the family of L-functions associated to Siegel modular forms, we assume the 
hp behave as independent random variables, each of which has a continuous probability 
distribution which is supported on [—10, 10] and which is bounded by 1, say. If we focus on 
one coefficient, say 63, and first choose all the other then fl4.11l) becomes 

C - 0.188349 63 = e X 0.00077, (4.12) 

where C is some number. Hence, there is a C so that, 

63 e [C - 0.004088, C + 0.004088]. (4.13) 

Since the PDF of 63 is assumed to be bounded by 1, the probability of (14.131) being true is 
less than the length of that interval, which is 0.00817. In other words, there is less than a 
1 percent chance that LpiQ would accidentally pass that test for satisfying the functional 
equation. We have proven: 

Theorem 4.2. Fix the parameters in the functional equation of L{s), as described at the 
beginning of Section 4-1 < o.iT'd suppose coefficients bj, j & J are known and the remaining 
coefficients obey the Ramanujan bound. Let L{s, (3) be given by (14. 8p . and choose real numbers 
Pi, P2 CLnd a complex number Sq. Write 

Z(so,/3i) -Z(so,/32) = ^Vjhj 

j 

= J2v,b, + e5, (4.14) 

where |B| < 1 and 6 is determined as described in Section 

If the bj forj G J are chosen independently from continuous probability distributions whose 
PDFs are bounded by 1, then the probability that (14.141) is consistent with the functional 
equation is less than 5/\vj\ for any j E J. 

It is easy to extend Theorem 14.21 to the case of several equations. Suppose we know bj for 
j G J, and choose Sk, Pk,i, and Pk,2 for 1 < k < K. We have 

Z{sk,/3k,i) - Z{sk,/3k,2) = Vkjbj 

j 

Select , . . . , , and suppose all the other bj have been determined. Then the system 

{Z{sk, /3k,i) - Z{sk, Pk,2) = 0}f=i (4.16) 

is equivalent to 

A{hj,,. . . , hj^) e (Ci, ...,Ck) + [Si, 5i] X • ■ ■ X [-6k, Sk], (4.17) 
where A is the matrix {vkjf.)- So we can rewrite the condition on the bj^. as 

{b,„. . . , h,,^) eA-\Ci,..., Ck) + A-\[-5i, 5i] X ■ ■ ■ X [-5k, 5k]). (4.18) 
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Since the PDFs of are assumed to be bounded by 1, the probabihty that f l4.18p occurs is 
bounded by the volume of the right side, which is 2^| det ■ ■ ■ 6k- 

Here is an example using some of the data from Table l4?n Pairing the 1st and 5th entries 
of column 2, and the 5th and 6th entries in column 3, we find: 

Z{\ + z, 0) - Z{\ + z, 1) = - 0.07393 + 0.0586962 + O.IOI7563 + ■ ■ ■ ± 0.00013 
Z{\ + 2i, 1) - Z{\ + 22, 2) = - 0.41376 + O.I8O2I62 + 0.4340163 + ■ ■ ■ ± 0.0077. (4.19) 

Using the coefficients of 62 and 63 in fl4.18p we find that the probability of the above system 
being satisfied for random 62? 63 is less than 0.00059. This strikes us as rather convincing 
evidence that the expected functional equation of Lp^iQ is in fact correct. 

We have a few comments on these calculations. Our purpose is to show that it is possible 
to quantify the precision to which changes to the test function in the approximate functional 
equation give a check on the functional equation of an L-function. Since it is, in a sense, 
nonsensical to treat the known coefficients of an L-function as random, we have tried not to 
push the analogy too far. If the coefficients hj really were independent and random, then the 
sum involving hj in equations like (14. lip would have a very large variance and the probability 
that (14. lip holds would be much smaller than our estimate. We chose to focus on just one 
or two coefficients at a time in order to not stretch plausibility too much. Also, the method 
of equation ( I4.18P gives poor results if the matrix is close to singular. In fact, one can add 
a new equation and increase the probability that the system is consistent, which is absurd. 
This can happen in practice: adding another equation based on Table 14.11 to ( I4.19P actually 
gives worse results. This is due to a dependence among the equations, arising from the fact 
that for small t it takes relatively few coefficients to evaluate L(| + zt). We will return to 
this topic in a subsequent paper. 

4.3. Rigorously evaluating L-functions. In Section HTT] we estimated the contribution of 
the terms involving the 6„ which were not known explicitly, but only assumed to satisfy the 
Ramanujan bound. This involves estimating the contribution of infinitely many terms and 
occurs in two steps. First, using Lemma Hl3] below, we determine and 61 so that the terms 
involving 6„ with n > N contribute, in total, less than 61. Then we explicitly evaluate the 
contributions of the terms fi{s,n) and /2(s,n) occurring in (14. 3 P for 83 < n < A^; call that 
contribution 62- Then our estimate for the contribution of the unknown terms is 61 + 62- 

For example, let /3 = 1 and s = | + i, in order to obtain the entry in the 5th row and 
second column of Table HH With A^ = 10, 000 we find 61 < 10"^ and 62 < 0.000373, as 
reported. This approach was used to determine the values in Table 14.11 and elsewhere in 
Section HI 

The following is a very slight modification of Lemma 5.2 of Booker 
Lemma 4.3. Let 

G*iu;vAH}) ■■=7^ [ e(^+^^^)(^-^)nrM(. + /x,)-. (4.20) 

Then for X > r, 

i=i 
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where6=^{l-\r]\), = \{^^^ - I) , ^ = \ + l{l + Y. l^j) , K = 2 



2r+l g(5(r-l) _TTrrj^j£ 



e 4 , and 



Note that our G* is identical to the function G in Lemma 5.2 of [1] except for the extra 
factor of 1/s in the integrand. 

Proof. Move the hue of integration to the 2cr hne and let s = 2a + 2it and use the trivial 
estimate l/|s| < l/{2a) to get 



G*(n;r?,{/i,})<^i- / 

271 J 2a 



dt. 



(4.22) 



Now exactly follow the proof of Lemma 5.2 in [1], which in the last step chooses a = X/r. □ 
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